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The existence of a second bound state of PsH that is electronically stable and also stable against 
positron annihilation by the normal 27 and 87 processes is demonstrated by explicit calculation. 
The state can be found in the ■^''*S° symmetries with the two electrons in a spin triplet state. The 
binding energy against dissociation into the H(2p) + Ps(2p) channel was 6.06 X lO"** Hartree. The 
dominant decay mode of the states will be radiative decay into a configuration that autoionizes 
or undergoes positron annihilation. The NaPs system of the same symmetry is also electronically 
stable with a binding energy of 1.553 x 10"'^ Hartree with respect to the Na(3p) + Ps(2p) channel. 

PACS numbers: 36.10.-k, 36.10.Dr, 34.85.+X 



The stability of a bound state composed of two elec- 
trons and a positron, the positronium negative ion, was 
first demonstrated in a seminal calculation by Wheeler 
[ij. Shortly after this calculation, the four body systems, 
PsH and Ps2 were shown to be stable 0, [3| ■ Since that 
time, only a few other electronically stable states have 
been discovered that can be formed from combinations 
of p\f-~ and e+. These are additional bound states of 
a compound that is best described as e+PsH 
0, and a 4e~, 26+) complex 0- Additionally, a 
number of atoms have been identified as being capable 
of binding positronium and positrons 0, [lo| 

A common feature of all these systems is that the 
positron annihilation process occurs by either a 27 or 87 
process with rates of order 10^ s~^ or 10^ s~^ (for those 
systems for which an annihilation rate has been deter- 
mined). In the present letter, we identify a new class of 
positronic compounds that are electronically stable, and 
in addition they have the unusual feature of decaying 
very slowly by 27 or 87 annihilation. Stable variants of 
PsH, and NaPs are identified and initial estimates of their 
binding energies are given. The existence of a new bound 
state of PsH is somewhat surprising given the amount 
of activity involved in ideritifying the resonant states of 
the Ps-H complex [Til H iH- The new PsH and NaPs 
bound states are unnatural parity states with symmetry 
conditions that act to prevent positron annihilation and 
to also prevent decay into the the lowest energy dissoci- 
ation products. These systems have the two electrons in 
a spin triplet state, a total orbital angular momentum of 
zero, and an odd parity, i.e. — 0~. Positron anni- 
hilation by the 27 or 87 process is forbidden for such a 
state. 

First consider the 27 process (which occurs at a rate of 
8 X 10^ s~^ for the Ps ground state). For this process to 
occur, the annihilating electron-positron pair must be in 
a spin singlet state and the relative angular momentum 
must be zero. (The decay rate is not absolutely zero since 



the Ps(2p) levels can undergo 27 and 87 annihilation at 
rate proportional to and respectively ^14, 15]. The 
rates for the different Ps(2p) levels have been calculated 



to be approximately 10 s 



115|.) 



Now consider the electron-positron annihilation of a 
PsH state of ^S° symmetry. The relative angular mo- 
mentum of the annihilating pair (iid) must be zero. This 
means the total angular momentum of the state will come 
from the center-of-mass motion of the annihilating pair 
(Lcm)j and from the angular momentum of the spectator 
electron (^spectator)- The total parity of the state is de- 
termined by the parity of the individual constituents, i.e. 

H = ( — 1) ^spectator +icm+irol^ jg glmply UOt pOSSlblC tO 

form an odd parity state with a total angular momentum 
of zero if any one of the constituent angular momenta is 
zero. Consequently, a two electron/one positron state of 
^S° symmetry cannot decay by the fast 27 process. 

These arguments also apply to the 87 annihilation pro- 
cess. The 87 process occurs for electron-positron pairs in 
a spin triplet state with a relative angular momentum of 
zero. Once again, it is simply impossible to form a state 
of ^S° (or "^S") symmetry if the relative angular momen- 
tum of the annihilating pair is zero. So it is reasonable 
to conclude that the lowest order 87 decay is not possible 
from a ^'"'S" state. 

These conditions also act to prevent the dissocia- 
tion of these four-body systems into combinations of the 
lower energy dissociation channels. Once again consider 
a ^S° state of PsH. Dissociation into Ps(ls)-|-H(ls) is 
forbidden since H = (—1)''^ where L is the orbital an- 
gular momentum between the Ps(ls) and H(ls) frag- 
ments. Similarly, dissociation into Ps(ns)-|-H(ni?) or 
Vs{ni)-\-Y{(ns) does not occur since it is not possible to 
construct an = 0^ state if one of the angular mo- 
mentum is zero. The lowest energy dissociation channel 
would be into Ps(2p)-|-H(2p) (p-wave) with an energy of 
—0.1875 Hartree. Another possible decay would be into 
the H~(2p^ '^P'^)+e+ channel with a threshold energy 



of -0.125355 Hartree [H, [B, [3- The stability of the 
H~(2p^ ■^P°) bound state also suggests a mechanism for 
binding. One can think of the positron trapped into a 2p 
state of the H~ attractive potential well. If the H~ state 
is regarded as a point particle with an internal energy of 
« —0.125 Hartree, then a positron in the 2p state will 
lower the total energy to —0.250 Hartree. In actuality 
the H~(2p^ ^P'') state is very diffuse, but this model does 
suggest that there is a large energy advantage associated 
with binding the positron to the negative ion. 

All the calculations in the present paper were per- 
formed with a configuration interaction approach 1^, 2^, 
2l| . The CI basis was constructed by letting the two 
electrons (particles 1 and 2) and the positron (particle 
0) form all the possible total angular momentum Lt — 
configurations, with the two electrons in a spin-triplet 
state, subject to the selection rules. 



max(£o,4,^2) < J 
in(^i,£2) 

(_l)('?0+«l+«2) 



min(£i,£2) < Ant , 



1 



(1) 
(2) 
(3) 



In these rules £o, and £2 are respectively the orbital 
angular momenta of the positron and the two electrons. 
We define {E)j to be the energy of the calculation with 
a maximum orbital angular momentum of J. The sin- 
gle particle orbitals were Laguerre Type Orbitals (LTOs) 
with a common exponent chosen for all the orbitals of 
a common e ^m,^- The orbitals basis sets for the 
positron and electrons were identical. 



TABLE I: The energy of the ^•''8° state of PsH as a function 
of J. The threshold for binding is —0.1875 Hartree. The 
column n gives the total number of occupied electron orbitals 
(the number of positron orbitals was the name) while Nci 
gives the total number of configurations. The results of the 
J — > 00 energy extrapolations at J = 10 are also given. 



J 


n 


Nci 


{E)j 


1 


15 


1800 


-0.16755817 


2 


30 


6975 


-0.17938456 


3 


45 


19125 


-0.18327387 


4 


60 


36000 


-0.18510508 


5 


75 


54675 


-0.18612672 


6 


90 


74925 


-0.18675211 


7 


105 


95175 


-0.18715811 


8 


120 


115425 


-0.18743280 


9 


135 


135675 


-0.18762315 


10 


150 


155925 


-0.18775631 


1- term eq.Q 

2- term eq.® 

3- term eq.® 






-0.18797567 
-0.18806917 
-0.18810659 



The Hamiltonian was diagonalized in a basis con- 



structed from a large number of single particle orbitals, 
including orbitals up to £ = 10. There were 15 radial 
basis functions for each £. Note, the symmetry of the 
state prevented the electrons or positrons from occupy- 
ing £ — orbitals. The largest calculation was performed 
with J = 10 and Lmt — 3 and gave a CI basis dimen- 
sion of 155925. The parameter Lint does not have to 
be particularly large since it is mainly concerned with 
electron-electron correlations The resulting Hamil- 
tonian matrix was diagonalized with the Davidson algo- 
rithm [22| , and a total of 300 iterations were required for 
the largest calculation. 

The energy of the PsH ^'^S° state as a function of J 
is given in Table ID The calculations only give an en- 
ergy lower than the H(2p) + Ps(2p) threshold of —0.1875 
Hartree for J > 9. A major technical problem af- 
flicting CI calculations of positron-atom interactions is 
the slow convergence of the energy with J [lol . 21 1. 
The J ^ 00 energy, {E)oo, is determined by the use 
of an asymptotic analysis. The successive increments, 
Ai?,7 = {E)j — {E)j-i, to the energy can written as an 
inverse power series [H, 2^ 24, 25, 2a|, viz 



AEj 



Bp 



Ce 



ij+^r (j+^y (j+^r 



(4) 



The first term in the series starts with a power of 6 since 
all the possible couplings of any two of the particles result 
in unnatural parity states 27 1. 
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FIG. 1: The binding energy, e = -((£) + 0.1875), of the 
^'*S° state of PsH as a function of J. The directly calculated 
energy is shown as the solid line while the J oo limits 
using eq. Q with 1, 2 or 3 terms are shown as the dashed 
lines. The H(2p) -I- Ps(2p) dissociation threshold is shown as 
the horizontal solid line. 



The J ^ 00 limit, has been determined by fitting sets 
of {E),j values to asymptotic series with either 1, 2 or 
3 terms. The coefficients, A^, Be and Ce for the 3- 
term expansion are determined at a particular J from 4 
successive energies {{E)j_3, {E)j_2, {E)j-i and {E)j). 
Once the coefficients have been determined it is easy to 
sum the series to oo and obtain the variational limit. 
Application of asymptotic series analysis to helium has 
resulted in CI calculations reproducing t he g round state 
energy to an accuracy of fslO"* Hartree [H, Si- 
Figure [1] shows the estimates of (£')oo as a function of 
J. A quick visual examination suggests that the extrap- 
olations are converging to a common energy while the 
energy of the three-term extrapolation is close to stabi- 
lized at J = 10. The impact of the extrapolations is sig- 
nificant since they more than double the binding energy. 
The most precise estimate of the binding energy is the 
three-term extrapolation at J = 10, namely 6.06 x 10~* 
Hartree (this is computed using the {E)oo in Table [T|. 



rect part of the core potential is attractive for electrons 
and repulsive for the positron. The impact of the direct 
and exchange part of the HF core interactions on the ac- 
tive particles are computed without approximation. One- 
and two-body semi-empirical polarization potentials are 
then added to the potential. The adjustable parameters 
of the core-polarization potential are defined by reference 
to the spectrum of neutral atom [i^, HO] ■ 

The system that was investigated was the ^'*S° state of 
NaPs. The energies of the 3s and 3p states in the model 
potential were -0.18885491 and -0.11156294 Hartree. 
The experimental binding energies are —0.188858 and 
—0.111547 Hartree respectively [sij. Electronic sta- 
bility requires a total 3-body energy of —0.17405849 
Hartree. The energy of the "^P*^ excited state of Na~ 
is —0.11342529 Hartree, i.e the Na(3p) has an electron 
affinity of 0.001862 Hartree with respect to attaching an 
electron to the ^P° state. This is reasonably close to the 
original value of Norcross, 0.00228 Hartree (29l |. 



TABLE II: The energy of the ^'"8° state of NaPs as a func- 
tion of J. The threshold for binding is -0.17410932 Hartree. 
The column n_ gives the total number of occupied electron 
orbitals, n+ gives the number of positron orbitals, while Nci 
gives the total number of configurations. The results of the 
J —> oo energy extrapolations at J = 10 are also given. 



2.0 



J 


n_ 


n+ 


Nci 




{E)j 




1 


16 


15 


2040 


-0 


153740547 


-0.02031794 


2 


31 


30 


7440 


-0 


166301513 


-0.00775597 


3 


46 


45 


19815 


-0 


170503790 


-0.00355470 


4 


61 


60 


36915 


-0 


172480804 


-0.00157768 


5 


76 


75 


55815 


-0 


173572387 


-0.00048610 


6 


91 


90 


76290 


-0 


174230717 


0.00017223 


7 


106 


105 


96765 


-0 


174650110 


0.00059162 


8 


121 


120 


117240 


-0 


174928718 


0.00087023 


9 


136 


135 


137715 


-0 


175119526 


0.00106104 


10 


151 


150 


158190 


-0 


175253262 


0.00119477 









1-term eq.Q 


-0.17547356 


0.00141507 


2-term eq.® 


-0.17556808 


0.00150960 


3-term eq.® 


-0.17561102 


0.00155252 



Having the established the stability of the ■^■''8° state 
of PsH, it is natural to ask whether other systems with 
this symmetry are stable. The obvious candidates are 
the alkali atoms, since some of them have np^ ^P® nega- 
tive ion bound states [1^ that can act as a parent state 
to bind the positron. The treatment of such systems re- 
quires the use of a frozen core approximation. The details 
of this approximation have been discussed in great detail 
elsewhere [il0,|2l|, so only the briefest description is 
given here. The model Hamiltonian is initially based on 
a Hartree-Fock (HF) wave function for the neutral atom 
ground state. The core orbitals are then frozen. The di- 
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FIG. 2: The binding energy (in units of Hartree) of the ^'"S" 
state of NaPs as a function of J. The directly calculated 
binding energy is shown as the solid line while the J — » oo 
limits using eq. @ with 1, 2 or 3 terms are shown as the 
dashed lines. The Na(3p) 4- Ps(2p) dissociation threshold is 
shown as the horizontal solid line, lines. 

The calculations upon NaPs were very similar in scope 
and scale to those carried out upon PsH. About the only 
difference was that an extra £ = 1 orbital was added 
to the electron basis. Table |lT] gives the 3-body energy 
(relative to the Na+ core) as a function of J. The binding 
energy ej is defined as ej = -({E) + 0.17410932). The 
positron complex is more tightly bound than for PsH and 
becomes electronically stable when J > 5. 

Figure [2] shows the variation of £oo as a function of J. 



4 



Once again the two and three term extrapolations seem 
to be converging to a common energy which is some- 
what larger than the best explicit calculation. The 3- 
term value of Soo determined at J = 10 was 0.001553 
Hartree. This is probably the the best estimate of the 
binding energy of the complex. The positron can annihi- 
late with the core electrons via the 2j process since the 
symmetry considerations are irrelevant here. However, 
the annihilate rate of Fcorc — 1-66 x 10^ s^^ is small 
because the positron cannot occupy a ^ = orbital. 

The PsH and NaPs ^■''8° complexes are stable against 
auto-ionizaton, and only decay slowly by positron anni- 
hilation. However there are other possible decay modes. 
Both these complexes can emit a photon, decaying to a 
state of ^^^P'^ symmetry. For example, a Ps(np) fragment 
in the complex can emit a photon decaying to a Ps(ls) 
type fragment. The Ps(ls) fragment could then annihi- 
late by the 2j or 37 process. In addition, a 2:4pc gi^g^i^g 
could also decay by auto-ionization. Due to their low 
binding energies, these systems can be expected to have 
a structure composed of an Ps(2p) cluster loosely bound 
to an atomic X(np) excited state. The lifetime of these 
states can be expected to be comparable to the lifetime 
of the fragments against single photon decay, e.g. H(2p) 
H(ls). So the overall lifetimes of the states can be 
expected to be of order 10~® - 10~^ seconds. 

It is possible that there are other positronic complexes 
of ^'"^S" symmetry that are bound. The K~, Rb^ and 
Cs~ ions have all been predicted to have np^ bound 
states. So the existence of a stable ^''^S" positronic com- 
plex would seem to be highly likely. It would also be 
interesting to examine Li as this does not appear to have 
a bound ^P'' negative ion [2^. 

Besides the alkali atoms, another physical system pos- 
sibly admitting an unnatural parity bound state would 
be the di-positronium molecule. There have been two at- 
tempts to find such a bound state, they were unsuccess- 
ful or inconclusive 
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33l |. However, the investigation of 
Bao and Shi showed that a ^S° state was very close to 
being bound, even if it was not bound [s^. This raises 
the tantalizing possibility that a more exhaustive calcu- 
lation might reveal the existence of a Ps2 state that de- 
cayed very slowly by positron annihilation. Besides the 
Ps2 molecule itself, there is the possible existence of a 
new biexciton excited state 3J]. The parent Ps^ ''P 



ion 

+ /m^- mass ratios 



is known to be stable for certain m 
[35I . [3^. In circumstances where the mass ratios make 
the '^P® state of the charged exciton {e~ ,e~ ,h) state sta- 
ble, it could be expected that a biexciton state of i'3,5go 
symmetry would be electronically stable. 

These calculations were performed on Linux clusters 
hosted at the South Australian Partnership for Advanced 
Computing (SAPAC) and SDSU Computational Sciences 
Research Center, with technical support given by Grant 
Ward, Patrick Fitzhenry and Dr James Otto. The au- 
thors would like to thank Dr D M Schrader and Dr C W 
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